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Abstract — We demonstrate that channel estimates 
obtained using the conventional correlation-based ap- 
proach are biased when short spreading sequences are 
employed in CDMA systems. A simple decorr elating 
estimator is proposed to remove the bias and hence 
improve estimation performance. 

I. Introduction 

Channel estimation is required for coherent detection in most 
code division multiple access (CDMA) receivers, including the 
Rake detector [1], the blind adaptive detector of [2], and mul- 
tistage interference cancellers (e.g. (3j). Inaccurate channel 
estimation will adversely affect their performance, especially 
the multiuser detectors 1 [4, 5]. 

In practice, the channel response is usually estimated with 
the aid of pilot symbols, through correlation with a pilot 
spreading code synchronized to the multipath delays, which 
are assumed to be estimated using a separate delay-tracking 
device. More complicated methods, based for instance on sub- 
space estimation (e.g. [6]), have also been proposed but will 
not be considered. 

In this paper, we show that correlation-based channel esti- 
mates are biased when short spreading codes (i.e. those with 
a period of one symbol interval) are used- Short codes are be- 
coming important 1 because they make certain multiuser de- 
tection techniques less complex (e.g. the decorrelator) while 
others (e.g. those based on identifying signal and noise sub- 
spaces |8]) need the time-invariant or cyclo-stationarity prop- 
erty of short codes to function. While detectors implement- 
ed using adaptive niters [7] do not require channel estimates, 
most others do and therefore channel estimation bias is a po- 
tentially serious problem that should be studied more thor- 
oughly. 

To remove the bias, we introduce a simple decorrelating 
estimator that requires an additional L 2 multiplications com- 
pared to the conventional estimator, where L is the number 
of paths tracked by the receiver. Simulation results show that 
the proposed estimator reduces channel estimation error and 
improves bit error rate (BER) substantially. 

II. Time Multiplexed Pilots 

We consider the synchronous downlink channel in this section, 
but it should be apparent shortly that our results hold for the 



asynchronous uplink as well. The DS-CDMA signal (sampled 
at chip rate) transmitted on the kih code channel is given by 
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*k(n) = ^2 dk(i)s k {n — iN) 



(D 



where d*(i) is the fcth symbol stream, s*(n) is the short 
spreading code used in the fcth channel and is non-zero on- 
ly in the range 0 < n < N> N is the processing gain and M is 
the total number of symbols transmitted in a packet. There 
are K code channels in total, and without loss of generality, 
the first one s\ («) carries the pilot symbols, meaning that P 
of the M symbols di(i) are known to the receiver. The various 
code channels may be assigned to different users, or to only 
one high-rate user, without affecting our conclusions. 

The channel is assumed to have L paths, characterised by 
time delays ri (which are in terms of sampling intervals) and 
complex attenuations q(i), / = 1, . . . ,i. This representation 
does not restrict the applicability of correlation- based channel 
estimation to slowly fading channels. However it does high- 
light the fact that this approach to channel estimation can de- 
liver only one estimate per symbol interval - if the maximum 
Doppler frequency is comparable to the symbol rate, then we 
will effectively be estimating the average or aggregate channel 

coefficient in each symbol interval. So in the fast-fading case, 
the channel coefficient Ci(i) should be understood to represent 
an average value. 

Collecting all received signal samples in the vector r, we 
can write 

L M~t 

r = £I><(') x <W + n (2) 

l = \ i=0 

where n is a zero- mean Gaussian noise vector with covariance 
matrix NqI, 



x / (i) = ^d>(i) s k =H 

jfc=l L 0( Af i)yv-r, J Jb=l 



<f*(i)s*,(0, (3) 



St = (sfc(0),6fc(l), . . . , Sk(N — 1)], Op is the p-dimensional vec- 
tor of zeros and s fc <(a) is implicitly defined in (3). 

The conventional method of channel estimation proceeds in 
two steps. First, rough estimates axe derived through correla- 
tion with the spreading code of the pilot channel, synchronized 
to the time delay of the desired path. Then, these rough es- 
timates are passed through a smoothing filter (for instance a 
moving-average filter) to produce smoothed estimates. 

The rough estimate of the first channel coefficient ci (i) is 



ci«) -dl(a)sri(t)r 



(4) 
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where d\(i) denotes the complex conjugate of di(i), a pilot 
symbol. Substituting (3) into (2), and ignoring inter-symbol 
interference terms, yields 

L 

s?,(t)r = 5Zci(i)di(i)s 1 , 1 (i)s„(i) 
/=i 

L K 

+ 53 X^ c,( ^ dfc W s " («)s*i(t) + *(z) (5) 



/=1 A =2 



w here = s^(t)n. Therefore 



c»(i) = C!(t) + y^c»(«)Pii(^0 

L K 

+ X^ZI^W^W^Wpi^i.O + rflW^t), (6) 

i = l it —2 

where p ife (m,n) = s** m (i)sfc»(0, the cross-correlation be- 
tween the mth and nth paths of the jth and A*h codes. We 
also assume without loss of generality that ||sn(*)|| = 1, and 
note that all spreading- code correlations are time- invariant. 
From (6), it is clear that 

L 

E[ci{i)\ = Cl (i) + 5^c,(0pn(l,0 ^^(i) (7) 

r=2 

or in other words the rough channel estimates, and hence the 
smoothed ones derived from them, are biased. When each 
user transmits over an independent channel, as on the uplink, 
equation (7) still holds because all ISI and MAI terms have 
zero mean. Similarly, when ISI is accounted for, the result 
applies because ISI has zero mean. 

III. I/Q Code Multiplexed Pilot Symbols 

In some cases, pilot bits are transmitted on the quadriphase 
(Q) carrier, while data is transmitted on the in-phase (I) car- 
rier. Separation of the two channels is accomplished through 
the use of orthogonal spreading codes for the I and Q channel- 
s. Other orthogonal code channels may be used in parallel by 
the same user, and of course other users can be transmitting 
simultaneously on other channels. 

As MAJ terms all have zero mean, they do not contribute 
to any bias in the channel estimates and it is therefore suffi- 
cient for our purposes to consider a simple scenario in which 
only one user transmits through a multipath channel. The 
transmitted signal will then be 

M-l 

*(«) = MO*i(« - iN) + jb p (i)s P (n - iN) (8) 

i=0 

where 61 (i) and 6 P (t) denote the data and pilot bits respective- 
ly, and 5i (n) and s p (-n) are the data and pilot spreading codes 
respectively, which are orthogonal over one symbol interval. 
The received signal vector will be 

L M — l 

r = S IM«)»m(«) + iM»)s pl (t)J + n (9) 



where 
su(t) — 



si ands P i(f) = Sp 

L o (A ,_<_ 1)W _ Ti j [ 0 (M _,_, 



si and s p are defined similarly to s fc in the previous section. 
A correlation-based estimate of ci will be 



(10) 



because it may be shown that in the absence of noise and 
multipath interference, ci (1) calculated in this way will equal 
ci . If we ignore ISI (which will have a mean of zero and not 
contribute to estimate bias anyway), 

L 

cx(i) = c,(i) + ^c(i)p 1 ,p(l,0 

1=2 

L 

+ 53 -**»M0*i (i)PD (1,0 + *(t). (11) 

Here : p pp (l,0 = s^s^i) and p P ,(l,l) = (i)s u <i) 

The mean of the rough estimate ci (i) is therefore once again 
not equal to cj (i), and so is biased. 

Long Codes When long codes are used, the second term 
on the right-hand sides of (7) and (11) become zero-mean 
random variables, and hence the bias disappears. The channel 
estimate bias is therefore a phenomenon created by the use 
of short codes, with conventional correlation- based channel 
estimation. However the bias term increases estimation error, 
and hence improved performance can be obtained even in long- 
code systems using the bias removal method discussed next. 

IV. A Bias-Free Estimator 
From (7) and (11), we see that in both TDM and I/Q code 
multiplexed cases, the mean value of the conventional channel 
estimate is 

^[c(0] = Rpc(i) (12) 
where c(t) s [cj(t), c 2 (i), . . . , Ci(»)J T is the vector of conven- 
tional estimates, 



Ppp(2,1) 



Ppp(£>1) 



1 



p pp (L,L- 1) 



Ppp(2,X) 
1 



and c(t) = [ci (*),..., cz(t)] T . For the time-multiplexed pi- 
lot system of section II, Pp^za.n) should be replaced by 
pxi (m,n). 

To create an unbiased estimate from c(i), we can "decor- 
relate" the estimates to form 



a 0 (t) = R^cfi). 



(13) 



The mean value of co(») will be equal to c, and hence it is 
an unbiased channel estimate. For short codes, R,, is time- 
invariant and its inverse only needs to be computed once. 
Thus the additional complexity in performing the decorrelat- 
ing operation effectively comes only from the I? multiplica- 
tions needed in (13). 

For long codes, we will cither have to re-compute Rp at each 
symbol interval, or implement the decollation operation us- 
ing linear multi-stage successive interference cancellation [9]. 
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Figure 1: Bit error rates for Rake receivers with con- 
ventional and bias-free channel estimates, in a short-code 
system. The curve labelled "AWGN" is for an AWGN 
channel, without ISI. P is the number of pilot symbols in 
a frame- 



V. Simulations and Conclusions 

The exact bit error rate of a linear receiver in a time- invariant 
multipath AWGN channel is readily computed as the sum of a 
number of Q functions (see for instance [10, Chap. 3]). Using 
this technique, we obtained the BER curves of Rake receivers 
with conventional (biased) and bias-free channel estimators 
for the I/Q code-multiplexed pilot scheme by 

1. Averaging the channel estimates over P pilot symbols; 

2. Finding the BER using those channel estimates; 

3. Repeating the process over 50 independent frames, and 
computing the average BER obtained. 

For a system with a processing gain of 8, one user and 
a four-path channel having path delays 0, 2, 4 and 6 chip- 
s' the BER for P = 20 and P = 40 are shown in Figure 
1. The pilot spreading code was Walsh code 0, and the da- 
ta spreading code was Walsh code 1. A scrambling code of 
[— 1, +1,-4-1, — 1,-1, +1,4-1, —1] was applied to both data and 
pilot channels. In this case, we see a gap of over 2 dB between 
the performance using biased and bias-free estimates at a BER 
of 10~ 3 . Also, while the performance of the bias-free algorith- 
m improves with a larger number of pilot symbols per frame, 
the same cannot be said for the conventional algorithm, since 
the errors in the latter case are due mostly to estimation bias 
and not to noisy estimates. 

The difference in channel estimation bias is even more 
vividly brought out in Figure 2, where we compare the con- 
ventional and decorr elating estimators. 100 independent tri- 
als were conducted, each of which selected a different set of 
spreading codes and channel coefficients. A three-path chan- 
nel was used, and the processing gain was eight. The E^/Na 
was fixed at 10 dB throughout, and the error magnitudes were 
obtained by averaging over 500 pilot symbols in each trial. We 
note that the error in the conventional estimator case is some 
ten times larger than the decorrelating estimator's. 

A final note is that the performance differences between 
the two algorithms are not always as large as in Figure 1. 




Figure 2: Magnitude of the channel estimate errors for 
100 independent channels, using the conventional and 
decorrelating estimators. 



Depending on the channel coefficients and the pilot spreading 
codes used, the bias terms in (7) and (11) may be so small 
that no discernible differences in performance can be observed. 
However, the bias-free estimates always perform better than 
the conventional estimates. 

We have shown in this paper that when short spreading 
sequences are used in CDMA, channel estimation done in the 
conventional correlation- based way will lead to biased esti- 
mates. In simulations, we have observed that this bias can 
lead to a substantial performance penalty in a Rake receiver. 
The performance loss in a multiuser detector has not been 
quantified and may be investigated in the future. To remove 
this channel estimate bias, we proposed a simple scheme akin 
to decorrelating multiuser detection. FVom simulations, we 
saw that the new bias-free channel estimation algorithm re- 
duced the magnitude of estimation errors by a factor of about 
ten. 
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